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Experimental and theoretical results on chemical functionalization of graphene are reviewed. Us- 
ing hydrogenated graphene as a model system, general principles of the chemical functionalization 
are formulated and discussed. It is shown that, as a rule, 100% coverage of graphene by complex 
functional groups (in contrast with hydrogen and fluorine) is unreachable. A possible destruction 
of graphene nanoribbons by fluorine is considered. The functionalization of infinite graphene and 
graphene nanoribbons by oxygen and by hydrofluoric acid is simulated step by step. 
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I. INTRODUCTION 

Chemical functionalziation is one of the main methods 
to manipulate the physical and chemical properties of 
nanoobjectsii^i^i^i^i^ii^i^ii^ and to study mechanisms of 
interaction of the nanoobjects with their environment 
(e.g., to investigate the stability of desired properties 
with respect to oxidation, hydrogenation, etc.). Initially, 
graphene was used by theoreticians as a simple model to 
describe properties of nanotube o^^'^^i-'^^'^'^i^^'^^i"'^''' . After 
the discovery of graphenei^ and of its extraordinary elec- 
tronic propertiesi2ii2£i2iiS2 the chemical functionalization 
of graphene became a focus of especial interest in contem- 
porary chemical physics. The main motivations of these 
studies are (i) modification of electronic properties via 
opening the energy gap in the electron spectrum of single- 
and bilayer graphen o^^i^'^'^^i^^'^'^i^^i^^ ; (ii) potential use 
of graphene for hydrogen storag o^'^'^'^i'^^'^^'^'^i'^^'^^ ; (iii) 
decoration of various defects in graphen o^-'^'^^i'^'^'^^i'^^'^^ ; 
(iv) a way to make graphene magnetic for poten- 
tial use in spintronic o^^i'^^'^^i'^^'^'^i'*'^ ; (v) search for 
ways to produce "cheap graphene" by chemical re- 
duction of graphite oxide and manipulation of its 
electronio^^i^i^ and mechanical^i^i^ properties; 
(vi) the functionalization of graphene edges in graphene 
nanoribbon g^^-^^-^^-^'^-^^-^^-^''-^^-^^'^"-^-^-^^-^'^'^'^-^^-^^-^''-^^-^^ 
and their protection^; (vii) oxidation and cracking of 
graphene2£ as tools to create graphene nanostructures 
of a given shape^liZSiiI^. 

There are two main types of functionalization, namely, 
a covalent one, with the covalent bond formation, and 
non-covalent, due to the van der Waals forces only. Most 
of the works deal with the covalent functionalization, and 
only few of the m^^i^*''''^:'''^''^'':''^'''^ treat the non-covalent 
one which is, of course, not surprising. First, the cova- 
lent functionalziation results in much stronger modifica- 
tion of geometric and electronic structures of graphene. 
Second, the most of the density functional codes used 
now do not allow to include the effects of the van der 
Waals interactions^Sii^'^ playing a crucial role in the 
non-covalent functionalization. In those cases the lo- 
cal density approximation (LDA)^ is mostly used, in- 
stead of the standard generalized gradient approxima- 
tion (GGA}§- which gives usually rather good results for 



layer compounds, such as graphite, hexagonal boron ni- 
tride and M0S2, due to error cancellationSS.. 

Sometimes, the chemical functionalziation in graphene 
is related to the ionic chemical bond. There are 
the cases of graphene layers in graphite interca- 
lated compounds^i^ and graphene layers at metal 
surfaces2^>28>89>9i^^^^^^ffiaoo. This type of func- 
tionalization is important to study possible superconduc- 
tivity in graphene, by analogy with the superconducting 
intercalates CaCg and YbCgiSi J£2^ 

Here we restrict ourselves only by the case of cova- 
lent functionalization of graphene. We will start with 
the hydrogenation of graphene, as a prototype of chem- 
ical functionazliation (Section 3). In Section 4, other 
examples will be considered, including realistic models 
of the functionalization by diatomic molecules, such as 
O2 and HF. In Section 5 we will discuss effects of finite 
width on the functionalization of the graphene nanorib- 
bons (GNR). We will give a resume of the work and dis- 
cuss some perspectives in Section 5. 



II. COMPUTATIONAL METHOD 

Our calculations have been performed with the 
SIESTA codei23ii04j£5,i06 ^^^^^ generalized gra- 
dient approximation (GGA)*^ dFT and TrouUier- 
Martins^^^ pseudopotentials. We used energy mesh 
cutoff 400 Ry, and /c-point mesh in Monkhorst-Park 
scheme^°^. During the optimization, the electronic 
ground states was found self-consistently by using norm- 
conserving pseudopotentials to cores and a double-^ plus 
polarization basis of localized orbitals for carbon, oxygen 
and fluorine, and double-C one for hydrogen. Optimiza- 
tion of the bond lengths and total energies was performed 
with an accuracy 0.04 eV /A and 1 meV, respectively. 

In our study of geometric distortions induced by the 
chemisorption of single hydrogen atom which will be 
described below (Section 3) a supcrcell containing 128 
carbon atoms with fc-point 4x4x1 mesh is used. For 
the cases of maximal coverage of graphene by various 
chemical species (Section 4) the standard elementary cell 
of graphene with two carbon atoms is used, similar to 
the case of full hydrogenation^, with fc-point 20x20x2 
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mesh. To simulate interaction of bulk graphene with 
hydrofluoric acid the supercell with 8 carbon atoms is 
used, similar to our previous works^iiS^ii^ with A;-point 
11x11x1 mesh. To investigate chemical functionaliza- 
tion of graphene nanoribbons (Section 5) zigzag stripes 
with the width 22 and 66 carbon atoms are used, similar 
to Ref4fl, with 1x13x1 fc-mesh. 

For each step of the functionalization, its energy is de- 
fined as Echem,, = E„ - E„_i - Eads, whcrc E„ and E„_i 
are the total energies of the system for n-th and n-1- 
th steps and Eads is the total energy of the adsorbed 
molecule. In contrast with the standard definition of the 
chemisorption which we used in the previous works, this 
definition allows to evaluate energy favor or disfavor of 
each step separately. We choose Fiads as a total energy of 
the corresponding molecule in gaseous phase. This is the 
most natural choice to estimate stability of chemically 
modified graphene. This is the crucial issue, for exam- 
ple, for the case of graphite oxide where the desorption of 
hydroxy groups in the main obstacle to derive a "cheap 
graphene"—. Note that molecular oxygen exists in both 
triplet (spin-polarized) and singlet (non-spin-polarized) 
states, the former being more energetically favorable by 
0.98 eV (we have found for this energy difference the 
value 1.12 eV). We will use here, as well as in the previous 
work^ the energy of singlet O2; to obtain the data with 
respect to magnetic oxygen one needs merely to shift up 
the presented energies by 1.12 eV. When considering ox- 
idation of nanotubes both energy of single t 
and tripleljiiiiii^iiii^ oxygen are used. 



III. HYDROGEN ON GRAPHENE 

Hydrogenation of carbon nanoobjects is a sub- 
ject of special interest starting from the discovery of 
fuUerenesii^, due to potential relevance of this issue for 
the hydrogen storage problem and a general scientific 
importance for chemistry. Soon after the synthesis of 
the first fuUerenes theoreticalii^iii^ and experimentaliii 
works appeared on the minimal hydrogenated fuUerene, 
C6oH2- It was shown in these works that the most stable 
configuration of hydrogen on the fuUerene corresponds 
to the functionalization of the neighboring carbon atoms 
(1,2, according to chemical terminology, see Fig. [T|). Fur- 
ther investigations result in discovery of numerous sys- 
tems with larger contents of hydrogen, up to CeoHagiiS. 
Further hydrogenation of the fullerenes turns out to be 
impossible, due to their specific geometric structure. It 
is important that all known stable compounds CqqUx 
contain even numbers of the hydrogen atoms. Together 
with other known facts^iS it allowed to formulate the 
main principle for the hydrogenation of the fullerenes: 
hydrogen atoms can be bonded either with the neigh- 
boring or with the opposite atoms in carbon hexagons 
(1,2 or 1,4, in the chemical terminology). Another is- 
sue of special interest was magnetism of hydrogenated 
fullerenes caused by unpaired electrons on broken bonds. 




FIG. 1: (color online) Optimized geometric structure of 
C60H2 

This can be realized only for odd numbers of hydro- 
gen atoms per buckyball^^^. Unfortunately, these hydro- 
genated fullerenes are unstable and observed only as in- 
termediate products of some chemical reactions; for even 
numbers of the hydrogen atoms magnetism does not sur- 
vive and disappears with a time^iSSi. It is worth to note 
that the hydrogenated fullerenes themselves are quite sta- 
ble chemically and do not loose hydrogen even under the 
action of high pressures^^^. 

The discovery of single-wall carbon nanotubesi^ii^^ 
(further we will discuss only this kind of the nan- 
otubes) and successful attempts of their reversible 
hydrogenatio n^^^'^^^ , which makes the nanotubes 
more prospective materials for hydrogen storage than 
fullerenes, have inspired numerous theoretical models of 
the nanotube hydrogenation. Geometric structure of the 
nanotubes differs essentially from that of the fullerenes. 
Single-wall carbon nanotubes can be represented as a 
graphene scroll of length up to few microns and width 
1-1.5 nm. 

Importantly, armchair and zigzag nanotubes exist with 
different structure of the honeycomb net on their sur- 
face, and the nanotubes can be chiral which is essen- 
tial for their chemical functionalization''. To avoid these 
complications, graphene is frequently used to simulate 
chemical functionalization of the nanotube a^"'i'^''^'^i'^''^ . 
However, as was discussed in Rcf.^^, a curvature radius 
of typical nanotubes (1.0 to 1.5 nm) is comparable with 
the radius of geometric distortions induced by hydrogen. 
Effects of the curvature of nanotubes on the hydrogena- 
tion was studied in detail recentlyi2£ and appeared to 
be quite essential, the hydrogenation being different in 
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two directions, along and around the nanotubes. Influ- 
ence of different types of chirality on chemistry of nan- 
otubes and differences with graphene was discussed also 
in Ref^^. 

Let us discuss briefly the main results of these works. 
Potential barriers for the cases of adsorption of a sin- 
gle hydrogen atom on graphene and several types of 
nanotubes, as well as formation of hydrogen chains on 
graphene and around the nanotubes have been studied 
in Ref.— . Modification of the electronic structure of 
perfect graphene and graphene with the Stone- Wales de- 
fect at the adsorption of single hydrogen atom was dis- 
cussed in Refjii. The authors- investigated a relation 
of magnetism with chemisorption of hydrogen for perfect 
graphene and for two types of defects, that is, monova- 
cancy and interlayer carbon atom in graphite. Sluiter 
and Kawazoei^ used the cluster expansion algorithm to 
find the most stable configuration for complete coverage 
of graphene by hydrogen, to model the hydrogenation of 
nanotubes. This configuration corresponds to so called 
graphane which has been studied afterwards by the DFT 
calculations'^^'^^. Wessely et al^^ simulated the core level 
spectra of carbon for the case of chemisorption of single 
hydrogen atom on graphene, to interpret existing exper- 
imental data on hydrogenated nanotubes. 

Graphene was used also to model the hydrogenation 
of graphit o"'^^^'^'^^'^'^^'^'^^ . After the discovery of graphene 
several works have been carried out on the hydrogena- 
tion of graphene itself, to modify its physical properties. 
In Ref.— the chemisorption of single hydrogen atom was 
studied in more detail than before, as well as interactions 
of magnetic moments arising, due to unpaired electrons, 
at the hydrogenation of carbon atoms belonging to the 
same sublattice (it was shown that in this case the inter- 
action turns out to be ferromagnetic). At the same time, 
according to Ref.— the combination of three hydrogen 
atoms (two of them belonging to sublattice A and the 
third one to sublattice B) remains magnetic which may 
be important in light of discussions of potential mag- 
netism of carbon systems. 

Sofo and coworkers^^ considered theoretically a hypo- 
thetic material, graphane, that is, graphite with com- 
pletely hydrogenated carbon layers (in each layer, all hy- 
drogen atoms coupled with one sublattice are situated 
above, and with another sublattice - below the layer). 
This structure corresponds to weakly coupled diamond- 
like layers, with sp'^ hybridization, instead of sp^ in 
graphite. This change of hybridization results in opening 
a gap of order 3 eV in electron energy spectrum, due to 
a transformation of tt and tt* orbitals to a and a* (see 
Fig. [2]). The cohesive energy was found to be relatively 
small (of order 0.4 eV per hydrogen atom) which allows to 
make the process of hydrogenation reversible^^. Roman 
and coworkers^ii^ considered different configurations of 
hydrogen on graphene for the case of one-side hydrogena- 
tion and have found that configurations when hydrogen 
atoms are bonded with different sublattices (similar to 
Fig. [3J:) are the most stable. In our workSi some general 
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FIG. 2: (color online) Band structure of a single graphene 
layer. Solid red lines are g bands and dotted blue lines are tt 
bands. 



principles of hydrogenation of graphene have been for- 
mulated based on calculations for various configurations 
in a broad range of coverage. 

Carbon atoms in graphene are in the sp^ hybridization 
state when each carbon atom has three a and one tt or- 
bitals (Fig. [2]). In contrast with numerous polyaromatic 
hydrocarbons with localized single- and double bonds tt- 
orbitals in graphene are delocalized and all conjugated 
chemical bonds are equivalent. Chemisorption of hydro- 
gen atom means a break of one of the tt bonds and trans- 
formation from sp^ to sp'^ hybridization. At the same 
time, unpaired electron, one of those forming the 7r-bond, 
remains at the neighboring carbon atoms. This electron 
is smeared in one of the sublattices (Figs. [Us and [4^) 
and forms a magnetic moment. The distribution of spin 
density in a radius of 1 nm correlates with distortions of 
the crystal lattice (Fig. [4)d). However, this situation is 
not robust in a sense that the chemisorption of the next 
hydrogen atom bonds this unpaired electron and kills the 
magnetic moment?-. 

Thus, the first principle of the chemical functionaliza- 
tion is the absence of unpaired electrons, or, in other 
words, of dangling bands. As a consequence of this prin- 
ciple, chemisorption of functional groups on different sub- 
lattices are much more energetically favorable than on 
the same sublattice. 

The next principle is a minimization of geometric frus- 
trations. One can see in Fig. that the chemisorption 
of single hydrogen atom results in essential atomic dis- 
placements inside a circle of radius 5 A (two periods of 
graphene crystal lattice) and in smaller but still notice- 
able distortions with a characteristic radius 10 A (Fig. 
[5]). The chemisorption of single hydrogen atom leads to 
a strong shift up of the carbon atom bonded with the 
hydrogen and shifts down for two neighboring atoms. It 
is worth to notice that if one fixes positions of all car- 
bon atoms except nearest and next-nearest neighbors 
of the central atom^SS, it changes essentially the pic- 
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FIG. 3: A schematic representation of chemical bonds in (a) 
pure graphene (dashed lines correspond to delocalized dou- 
ble bands), (b) graphene with single chemisorbed hydrogen 
atom (white circle, the black circle corresponds to unpaired 
electron at the broken bond, it distributes in the sublattice 
shown in red), (c) graphene with chemisorbed pair of hydro- 
gen atoms (nonequivalent sublattices are shown in different 
colors), (color online) 
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FIG. 4: Charge distribution for unpaired electron as a func- 
tion of a distance from the central carbon atom where hy- 
drogen is chemisorbed (a) and deviations of carbon atoms 
from flat configuration for the cases of single hydrogen atoms 
(red solid line) and two atoms chemisorbed on neighboring 
sites from different side of the graphene sheet (green line) 
(b). (color online) 



ture of atomic displacements. Such a procedure may 
be reasonable to simulate chemical functionalization for 
constrained systems such as graphene on graphite^"*^, 
Rui2ii2^i21, Iri^ and P1*21 where part of graphene lattice 
is strongly coupled with the substrate. 

Chemisorption of the next hydrogen atom is more fa- 
vorable when it results in minimal additional distortions. 
This means that the most energetically stable configu- 
ration arises if two hydrogen atoms are chemisorbed by 
neighboring carbon atoms at different sides of graphene 
sheet. One can see in Fig. |4j3 that in this case no es- 
sential additional atomic shifts are necessary. If only one 
side is available for the chemisorption the most favorable 
configuration of two hydrogen atoms is the bonding with 
carbon atoms in the opposite corners of the hexagon in 
honeycomb lattice (cites 1,4, or para-position, accord- 
ing to the chemical terminology}^^. For the case of two 
more distant hydrogen atoms they should migrate to this 
optimal position overcoming some potential barriers^^ 
which was observed experimentally for both graphite^^^ 
and graphenei^. 

These two principles of the chemical functionalization 



predict that in the most stable configuration for hydrogen 
pair on graphene the hydrogen atoms will sit on neigh- 
boring carbon atoms at the opposite sides with respect to 
the graphene sheet. Due to geometric distortions a region 
of radius 5 A around the pair is more chemically active 
than pure graphene, similar to the case of graphene with 
defectai2ii2^i^i2Si^ and, thus, next hydrogen atoms will 
be chemisorbed near the pair. In the case of graphene 
with defects this process involves some potential barriers 
whereas for ideal infinite graphene it will proceed without 
any obstacles, up to complete hydrogenation and forma- 
tion of graphane. 

In further works'^'*' the results^^ on stability of non- 
magnetic pairs of hydrogen on graphene and on the for- 
mation energies have been confirmed, also the case of 
polyaromatic hydrocarbons C54Hig^ and C42Higi^^ and 
C96H24^'^^ being considered. Recent paper-i^i deahng 
with adsorption and desorption of biphenyl on graphene 
demonstrates importance of step by step simulation of 
the chemical reaction. 

It is instructive to compare peculiarities of chemi- 
cal functionalization of graphene and carbon nanotubes. 



FIG. 5: Regions of strong (blue circle) and weak (yellow cir- 
cle) distortions of graphene lattice at the chemisorption of 
single hydrogen atom (small red circle), (color online) 

Similar to graphene, in nanotubes the most favorable 
place for chemisorption of the second hydrogen atom be- 
longs to the opposite sublattice with respect to the first 
one^^^. There is no direct experimental evidences of this 
structure for the case of hydrogen, however, for another 
functional groups it is known that normally they are at- 
tached to pairs of neighboring carbon atoms^i^ii^. Min- 
imization of geometric frustrations plays an important 
role also in the case of nanotubes. In contrast with flat 
graphene, a particular shape of the nanotubes makes the 
chemisorption from external side more favorable, with 
the adsorbed pairs situated along the nanotube. For large 
enough groups they are usually attached to the opposite 
sides of the nanotube cross section^, again, to minimize 
the carbon net distortions. In contrast with graphene 
where complete one-side hydrogenation is unfavorable 
for the case of nanotubes it turns out to be possible, 
from external side, as is shown both theoretically^^ and 
experimentalljJ^ii2^. 



IV. FUNCTIONALIZATION OF GRAPHENE 
BY OTHER CHEMICAL SPECIES 

As discussed above, graphane is not very stable which 
is actually quite good from the point of view of hydrogen 



FIG. 6: Top and side view of graphene functionalized by CH3 
groups (optimized atomic positions). Carbon and hydrogen 
atoms are shown in green and white, respectively, (color on- 
line) 



storage since it allows to hydrogenate and dehydrogenate 
graphene at realistic temperatures^-. At the same time, 
for potential use of graphane in electronics this can be 
considered as a shortcoming. Therefore it is interest- 
ing to consider other functional groups to search a com- 
pound with the electronic structure similar to graphane 
but larger cohesive energy. 

The computational results are presented in Tableland 
in Fig. [HI Based on them one can divide all the functional 
groups in two classes, (a) with the energy gap of order 3 
eV and large enough cohesive energy and (b) with much 
weaker bonding and much smaller energy gap, or with- 
out a gap at all. This difference can be related with the 
length of chemical bonds between carbon atoms (carbon- 
carbon distance) . For the substances from the first class 
it is almost equal to that in diamond, 1.54 A whereas 
for the substances from the second class it is larger. This 
means that for the latter case, due to mutual repulsion of 
the functional groups, graphene is destabilized which di- 
minishes the cohesive energy. This situation is observed 
experimentally in graphite oxide which cannot be com- 
pletely functionalized due to interaction between hydroxy 



TABLE I: 

Dependence of the carbon-carbon bond length (d), in A , and 
the electron energy gap {AE), in eV, on chemical species for 
the case of 100% two-side coverage of graphene. 



Chemical species 


d 


AE 


NO 


1.622 


none 


NH2 


1.711 


2.76 


CN 


1.661 


2.11 


CCH 


1.722 


0.83 


CH3 


1.789 


0.84 


OH 


1.620 


1.25 


H 


1.526 


3.82 


F 


1.559 


4.17 


H and F 


1.538 


5.29 


H and CI 


1.544 


4.62 



groups^^, in agreement with the XPS datai^. 

As examples of chemical species which can provide, 
similar to hydrogen, a complete coverage of graphene, 
we will consider F2, HF and HCl. Fluorine seems to be 
very promising for the functionahzation of graphene since 
it should produce a very homogeneous structure, with a 
complete coverage by atoms of the same kind. Therefore, 
this might be a way to create a two-dimensional crystal 
with rather large energy gap but high electron mobil- 
ity, due to small degree of disorder. However, fluorine 
is very toxic and very aggressive which may be a prob- 
lem for industrial use. Also, fluorine is used for ripping 
of nanotube^i^ii^ which means that, potentially, edges 
and defects of graphene lattice interacting with the fluo- 
rine can be centers of its destruction. Actually, graphene 
samples are always ripple d^"^^'^'^^'^'^^'^'*^'^"^^'"'^^^ or can be 
deformed as observed for graphene ribbons2^. 

Let us consider now, step by step, interaction of 
graphene with inorganic acids HF and HCl. The process 
of functionahzation of graphene by HF starts with the 
chemisorption of fluorine and hydrogen atoms at neigh- 
boring carbon sites at the same side of graphene sheet 

(Fig. m- 

As a first step of chemisorption of hydrofluoric acid 
molecule we consider bonding of hydrogen and fluorine 
atoms with the neighboring carbon atoms at the same 
side of graphene sheet (Fig. [7^). In contrast with molec- 
ular fluorine where the chemisorption energy is negative 
(see above) the energy of this step turns out to be pos- 
itive, -1-1.46 eV, mainly, due to strong distortions of ini- 
tially flat graphene at the one-side functionahzation. As 
the next step, we consider chemisorption of one more 
H-F pair from another side of graphene (fig. [7)d). Sim- 
ilar to chemisorption of the second hydrogen atom (sec 
above) this step is energetically favorable, with the for- 
mation energy -1.85 eV. Reconstruction to the structure 
shown in Fig. [T]; diminishes the total energy by 0.33 eV 
which makes this a possible third step of the process un- 
der consideration. After that, the chemisorption energy 
of third HF molecule turns out to be much smaller than 
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FIG. 7: A sketch illustrating processes and energetics of dif- 
ferent reactions of graphene with hydrofluoric acid. White 
and red circles are hydrogen and fluorine atoms, respectively, 
all energies are in eV. (color online) 



for the first one but still positive, -1-0.27 eV (Fig. [7Ji). 
Starting from the fourth molecule, the chemisorption en- 
ergies are negative so further reaction is exothermal. As 
a result, complete coverage of graphene by hydrofluoric 
acid is energetically favorable but requires high enough 
energy for the first step. This means that, the most prob- 
ably, bulk graphene without defects will be stable enough 
with respect to reaction with HF, at least, at room tem- 
perature and normal pressure. It seems to be practi- 
cally important since hydrofluoric acid is used to sol- 
vate Si02 substrate when prepare freely hanged graphene 
membranei^. One can expect in this case a formation 
of a chemically modiflcd layer near edges with essentially 
difl'erent electronic structure (see below). 

Our calculations show that the first step of functional- 
ization of graphene by hydrochloric acid is essentially dif- 
ferent, due to larger interatomic distance (1.27 A for HCl 
pair, instead of 0.97 A for HF) which close to carbon- 
carbon distance in graphene (1.42 A). As a result, there 
is no break of bonds between hydrogen and chlorine but, 
instead, HCl molecule hangs over graphene layer bond- 
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monolayer: -1-0.23 eV; bilayer: -(-0.48 eV 

-1^ 




FIG. 8: Optimized geometric structure for step-by-step ad- 
sorption of oxygen molecule on graphene. (color online) 
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ing with it by only weak van-der-Waals bond and there 
is no chemical reaction. Thus, the hydrocloric acid can 
be safely used to clean graphene. 

Let us consider now, using oxygen molecule as an ex- 
ample, a chemisorption of diatomic molecule with double 
chemical bond. In this case, at the first step the double 
bond will transform into the single one, with formation 
of single bonds with neighboring carbon atoms (Fig. [5]). 
This reaction requires to overcome some energy barrier 
which is, however, rather low since flexibility of graphene 
allows to minimize the energy costs due to lattice dis- 
tortions. Bilayer graphene is much less flexible, due to 
interlayer coupling, which makes the chemisorption en- 
ergy higher, in agreement with the experimental dataii^. 

If we would calculate chemisorption energy with re- 
spect to not singlet but triplet state of O2 molecules it 
would increase all chemisorption energies by 1.12 eV (see 
Section 2). The resulting energies are so high that it 
seems in contradiction with the experimental fact that 
graphene is relatively easily oxidized already at temper- 
ature up to 200 °G^. Also, "unzipping" process at 
the oxidation of graphite observed experimentall y ''"1^^° 
will involve energetically unstable states. Graphene ox- 
ide turns out to be also unstable with respect to triplet 
oxygen^^ii^. Probably these discrepancies mean that 
formation of singlet oxygen plays some role in these pro- 
cesses. This issue requires further investigation. A simi- 
lar problem was discussed earlier for nanotubeaiS^. 

Oxidation energy for graphene turns out to be lower 
than for the case of nanotube a^°^'^^° . In the latter case, 
TT-orbitals are rotated one with respect to another— which 
makes formation of oxygen bridges between carbon atoms 
more difficult. 

Previous computational results for oxidation of 
graphene^^^ are qualitatively similar to ours. However, 
the energy of cyclo-addition there is higher than ours. 
There are two possible reasons for this difference, the use 
of LDA in Reff^^ instead of GGA in our calculation and 
the use of relatively small graphene "molecule" (with 64 
atoms) instead of infinite system with periodic conditions 
here. 



OH NHNH2 



c) 

25% 



0.03 



-1.06 



FIG. 9: Sketch of chemical reduction of graphene oxideW? for 
various degrees of its functionalization. All energies are in eV. 
(color online) 



Let us consider now another important issue, the sim- 
ulation of reduction of graphite oxide (GO). We have 
proposed earlier— models of this compound with differ- 
ent degrees of its reduction. In ref.^'"^^ a process of the 
reduction has been studied and the scheme of correspond- 
ing reactions has been proposed. Using our model and 
this scheme, we have simulated three possible ways of 
the chemical reduction of GO depending on degree of 
preliminary reduction and calculated the corresponding 
energies (Fig. [S]). To compare, we have calculated also 
energy costs of direct reduction (e.g., by heating). One 
can see from the Figure that for high degrees of cov- 
erage of GO (weakly reduced) energy costs of chemical 
reactions are smaller than those of the direct reduction 
whereas for the case of strongly reduced GO, vice versa, 
direct total reduction turns out to be more energetically 
favorable. Similar, other reactions of real^^^ and poten- 
tial functionalization of GO can be simulated which may 
be important to search ways of its complete reduction to 
pure graphene. 
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FUNCTIONALIZATION OF GRAPHENE 
NANORIBBONS 



■^-Jt^-flHl-M-^l-(HH)-flHl-ONa-^^^ 



■2,50 



FIG. 10: A scheme of step-by-step oxydation of GNR of width 
2.2 nm. AU energies are in eV (top panel) . Optimized geomet- 
ric structure for partialy oxidized GNR (lower panel), (color 
online) 



Polyaromatic hydrocarbons (PAH)^— are quite close 
to graphene nanoribbons (GNR), and their chemistry, 
in particular, hydrogenation of their edges, is rather 
well studied™'. Another nanoobjects which can be con- 
sidered as predecessors of GNR are nanographite and 
nanographite intercalates which were studied, in partic- 
ular, in connection with their magnetisnii^'i^. Model 
calculations of their magnetic properties are presented in 
RefsJ^ii^ whereas first-principle calculations, using a 
fragment of graphene sheet as a model of nanographite, 
have been performed in Ref.-^^. In that paper, function- 
alization of the graphene edges by single hydrogen atom 
has been considered, as well as suppression of magnetism 
by fluorination, in relation with the experimental data-5^. 

Starting from the first work by S. Louie and 
coworkeraSi predicting, on the base of the DFT cal- 
culations, half-metallic ferromagnetic state of zigzag 
graphene edges, this problem attracts a serious 
attention^i^. Later, the first results have been reconsid- 
ered by the same group using more accurate many-body 
GW method^. Chemical fucntionalization of GNR is 
currently one of the most popular subjects in graphene 
science. 

O. Hod et al. have considered effects of the shape 
of semiconducting GNR on their stability and electronic 
structure^, electronic structure of graphenic nanodots^ 
(similar calculations of magnetism and electronic struc- 
ture have been presented also in Rcf.^^), effects of 
geometry of graphene nanoshcets on their electronic 
structure^ and enhancement of half-metallicity due to 
partial oxidation of the zigzag edges^. Interesting ef- 
fects which can arise in graphene nanoflakes are dis- 
cussed in Ref.^^. Usually, it is supposed in these cal- 
culations that edges of nanoflakes are passivated by hy- 
drogen atom o^^'^'^i^^'"'^^'-''^^"'^'"'^^^ but sometimes pure edges 
are consideredi^ii^. Real structure of some graphene 
nanoflakes was discussed in Refj^^. In general, the state 
of real graphene edges is not well understood yet and this, 
very important, problem needs further investigations. 

In Refi^ effect of vacancies at the edges on the elec- 
tronic structure of GNR were studied. Influence of 
chemisorption of different chemical species on the elec- 
tronic structure of GNR was investigated in Refsi ^^'^^ . 
Ref.^^ considered the chemisorption of molecular oxygen 
on GNR and its effect on their magnetism. Various con- 
figurations of hydrogenated§§ and self-passivatectfS, GNR 
were studied; in both cases nonmagnetic state of the 
GNR turned out to be the most stable. Refs^SiS demon- 
strated suppression of magnetism at the zigzag edges by 
oxidation. Since the latter is a natural process at the pro- 
duction of GNR^^iiiii^ some special measures to protect 
the magnetism chemically are probably necessary. 

Thus, works on the functionalization of GNR can be 
divided into three groups. First, there are investigations 
of electronic structure and physical properties (especially. 
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magnetism) of the edges themselves. This group includes 
most of theoretical works. Second, there are works study- 
ing the functionalization of the nanoribbons between the 
edges. It seems that the effect of the edges on the elec- 
tronic structure of GNR is long-range enough and, as a 
consequence, chemistry of GNR can be essentially differ- 
ent from that of bulk graphene. Research in this direc- 
tion is just starting. Third, there are attempts to find 
optimal geometric and chemical structure of GNR tak- 
ing into account realistic technological processes of their 
formation. Preliminary, there is an impression that some 
structures considered theoretically before can be very dif- 
ferent from the real ones, so we are, indeed, only in the 
beginning of the way. On the other hand, a lot of exper- 
imental and theoretical information is available for re- 
lated compounds such as aryneS'iSl, polyhexacarbens^^ 
and polycyclic aromatic hydrocarbons^^ which may be 
also relevant for GNR. 

It was shown in our previous work^, using hydrogen as 
an example, that the edges of GNR are centers of chemi- 
cal activity and their functionalization can be just a first 
step to the functionalziation of the whole nanoribbon. 
Here we will consider step by step the oxydation process 
of GNR. The computational results are presented in Fig. 
[TOl Thus, oxidation of GNR is a chain of exothermal pro- 
cesses. This result may be relevant also for better under- 
standing of burning and combustion of carbon systems. 
These issues, as well as more detailed results concern- 
ing oxidation and unzipping''*^ of graphene and related 
compounds will be discussed elsewhere. 

Let us consider now a reaction of GNR with the hy- 
drofluoric acid. A scheme of this reaction for the nanorib- 
bon of width 2.2 nm is shown in Fig. Illb . and depen- 
dence of the energy as a function of coverage is presented 
in Fig. [TTb . In contrast with the case of ideal infinite 
graphene for narrow enough GNR two steps with posi- 
tive chemisorption energy. This energy cost is so small 
that the process of hydrofluorination of this nanoribbon 
can take place even at room temperature. The situation 
is essentially different for the nanoribbon of width 6.6 nm. 
The steps with positive chemisorption energy just disap- 
pear, due to larger distance between the regions of the 
reaction. Thus, the hydrofluorination of a broad enough 
GNR is easy and complete. 

One can see, that, similar to the case of 
hydrogenationiSi (see also the previous section) GNR 
are very chemically active which probably means that, 
at least, some of their applications will require an inert 
atmosphere. 



VI. CONCLUSIONS AND PERSPECTIVES 

In Section 3 we discussed, using hydrogenation as 
example, three main principles of chemical function- 
alization of graphene: (i) broken bonds are very un- 
favorable energetically and, therefore, magnetic states 
on graphene are usually very unstable; (ii) graphene is 



very flexible, and atomic distortions influence strongly 
on the chemisorption process; (iii) the most stable 
conflgurations correspond to 100% coverage for two- 
side functionalization^^ and 25% coverage for one-side 
functionalization^. Based on these principles one 
can model the functionalization of not only perfect 
graphene2ii2^ii^ but also graphene with intrinsic and ex- 
trinsic defects which are centers of chemical activity""^. 

In Section 4 we have studied chemical functionaliza- 
tion of graphene by fluorine and hydrofluoric acids. They 
can provide complete coverage and, thus, semiconducting 
state with large enough electron energy gap and weak dis- 
order. Alternatively, non-covalent functionalization can 
be used to create the gap. 

The current situation seems to be a bit controversial. 
In some papers^^'^^ an energy gap opening due to ph- 
ysisorption of various molecules has been found, in con- 
trast with other resultsAS. i^'^/^i'^^/'^ . It is difficult to com- 
pare these works directly since they are, in general, done 
for different chemical species and their concentrations. 
However, keeping in mind that in some cases^S^^^i^l the 
same case of water has been studied one can conclude 
that the electronic structure is very sensitive to specific 
concentration and specific geometric configuration of wa- 
ter on graphene surface. At the same time, one should 
keep in mind that a very specific structure of water^i^S 
and other hydrogen-bonded substances may be not very 
accurately reproduced in the DFT calculations. Prob- 
ably, only combination of the DFT, quantum chemical 
calculations and molecular dynamics can clarify the sit- 
uation. 

Anyway, since the cohesive energy at physisorption 
does not exceed 20 kJ/mol the physisorbed graphene can- 
not be very robust which may result in some restrictions 
of its use in electronics. 

We have demonstrated using HF as an example that it 
is important to simulate the chemisorption process step 
by step since it allows to estimate accurately chemisorp- 
tion energies relevant for each step of the process. It turns 
out that for perfect infinite graphene the energy cost of 
first step for the case of HF is rather high whereas for HCl 
this reaction is practically impossible. This means that 
these acids can be safely used for cleaning of graphene 
samples. At the same time, the potential barrier for the 
oxidation is rather low which makes even moderate an- 
nealing of graphene in oxygen-contained atmosphere not 
very safe procedure. 

We have shown also, in the last section, that the 
graphene nanoribbons are chemically active enough to 
be oxidized and to interact with hydrofluoric acid with 
relatively low activation energy. It means that to prepare 
the graphene nanoribbons with a given chemical compo- 
sition alternative ways should be used, e.g., its synthesis 
from polyaromatic hydrocarbon o^^^i^^^ or even use mi- 
croorganisms as it is used already for biogenic growi22ii2i 
and cleaningii^. 
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FIG. 11; A scheme of step-by-step reaction of GNR of width 
2.2 nm with HF; the energies are in eV (a). The same energies 
as functions of coverage for the case of GNR with width of 
2.2 nm (red solid line) and 6.6 nm (green dashed line), (color 
online) 



